Adiabatic-Impulse approximation for avoided level crossings: from phase transition 
dynamics to Landau-Zener evolutions and back again 
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We show that a simple approximation based on concepts underlying the Kibble- Zurek theory of 
second order phase transition dynamics can be used to treat avoided level crossing problems. The 
approach discussed in this paper provides an intuitive insight into quantum dynamics of two level 
systems, and may serve as a link between the theory of dynamics of classical and quantum phase 
transitions. To illustrate these ideas we analyze dynamics of a paramagnet-ferromagnet quantum 
phase transition in the Ising model. We also present exact unpublished solutions of the Landau-Zcner 
like problems. 
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I. INTRODUCTION 

Two level quantum systems that undergo avoided level 
crossing play an important role in physics. Often they 
provide not only a qualitative description of system prop- 
erties but also a quantitative one. The possibility for 
a successful two level approximation arises in different 
physical systems, e.g., a single one-half spin; an atom 
placed in a resonant light; smallest quantum magnets; 
etc. 

In this paper we focus on avoided level crossing dynam- 
ics in two level quantum systems. Two interesting related 
extensions of above mentioned long list are the critical 
dynamics of quantum phase transitions P, S S S 
and adiabatic quantum computing . In both of these 
cases (and, indeed, in all of the other applications of the 
avoided level crossing scenario) there are many more than 
two levels, but the essence of the problem can be still cap- 
tured by the two level, Landau-Zener type, calculation. 

Our interest in Landau-Zener model originates also 
from its well known numerous applications to different 
physical systems. In many cases there is a possibility 
of more general dependence of the relevant parameters 
(i.e. gap between the two levels on time) than in the 
original Landau-Zener treatment. This motivates our ex- 
tensions in this paper of the level crossing dynamics to 
asymmetric level crossings and various power-law depen- 
dences. Appropriate variations of external parameters 
driving Landau-Zener transition can allow for an experi- 
mental realization of these generalized Landau-Zener like 
models in quantum magnets and optical lattices @ . 

We focus on evolutions that include both adiabatic and 
diabatic (impulse) regime during a single sweep of a sys- 
tem parameter. For simplicity, we call these evolutions 
diabatic since we assume that they include a period of 
fast change of a system parameter. This is in contrast to 
adiabatic time evolutions induced by very slow parameter 
changes when the system never leaves adiabatic regime. 
We will use a formalism proposed recently by one of us 
[3- It originates from the so-called Kibble-Zurek (KZ) 
theory of topological defect production in the course of 
classical phase transitions [iflllji ^^'^ works best for two 



level systems. The two level approximation was recently 
shown by Zurek et al and Dziarmaga ^ to be useful in 
studies of quantum phase transitions confirming earlier 
expectations - see Chap. 1.1 of Q and Q. Therefore, we 
expect that the formalism presented here will provide a 
link between dynamical studies of classical and quantum 
phase transitions. 

The present contribution extends the ideas presented 
in 0. In particular, we have succeeded in replacing the 
fit to numerics used there for getting a free parameter 
of the theory, with a simple analytic calculation. We 
also show that the method of 9] can be successfully ap- 
plied to a large class of two level systems. Moreover, we 
present nontrivial exact analytic solutions for Landau- 
Zener model and the relevant adiabatic-impulse approx- 
imations not published to date. 

Section|n]presents basics of what we call the adiabatic- 
impulse approximation. Section IIIII shows adiabatic- 
impulse, diabatic and finally exact solutions to different 
versions of the Landau-Zener problem. In Section Hvl we 
discuss solutions for a whole class of two level systems 
on the basis of adiabatic-impulse and diabatic schemes. 
Section Ivl presents how results of Section UTTl can be used 
to study quantum Ising model dynamics. Details of ana- 
lytic calculations are in Appendix 1X1 fdiabatic solutions). 
Appendix IbI f exact solutions of Landau-Zener problems) 
and Appendix [CI (adiabatic-impulse solution from Sec. 
CVl). 



II. ADIABATIC-IMPULSE APPROXIMATION 

We will study dynamics of quantum systems by as- 
suming that it includes adiabatic (no population transfer 
between instantaneous energy eigenstates), and impulse 
(no changes in the wave function except for an overall 
phase factor) stages only. The nature of this approach 
suggests the name adiabatic-impulse (AI) approximation. 
The AI approach originates from the Kibble-Zurek (KZ) 
theory of nonequilibrium classical phase transitions and 
we refer the reader to for a detailed discussion of AI- 
KZ connections. Below, we summarize basic ideas of the 
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KZ theory, and the relevant assumptions used later on. 

Second order phase transitions and avoided level cross- 
ings share one key distinguishing characteristic: in both 
cases sufficiently near the critical point (defined by the 
relaxation time or by the inverse of the size of the gap 
between the two levels) system "reflexes" become very 
bad. In phase transitions this is known as "critical slow- 
ing down" . 

In the treatment of the second order phase transitions 
this leads to the behavior where the state of the sys- 
tem can initially - in the adiabatic region far away from 
the transition - adjust to the change of the relevant pa- 
rameter that induces the transition, but sufficiently close 
to the critical point its reflexes become too slow for it 
to react at all. As a consequence, a sequence of three 
regimes (adiabatic, impulse near the critical point, and 
adiabatic again on the other side of the transition) can 
be anticipated. In the second order phase transitions this 
allows one to employ a configuration dominated by the 
pre-transition fluctuations to calculate salient features of 
the post-transition state of the order parameter (e.g., the 
density of topological defects). Key predictions of this 
KZ mechanism have been by now verified in numerical 
simulations |l2l ll3[ and more importantly in the labora- 
tory experiments |14L ITEl IT^ IT^ . 

The crux of this story is of course the moments when 
the transition from adiabatic to impulse and back to adia- 
batic behavior takes place. Assuming that the transition 
point is crossed at time t = 0, this must happen around 
the instants ±i, where t equals system relaxation time r 
As the relaxation time depends on system param- 
eter £, which is changed as a function of time, one can 
write 

rieii))=l (1) 

This basic equation proposed in Tl'l can be solved when 
dependence of r on some measure of the distance from the 
critical point (e.g., relative temperature or coupling e) is 
known. Assuming, e.g., r = tq/e and e = t/rq, where 
tq , the "quench time" , contains information about how 
fast the system is driven through the transition, one ar- 
rives at To/{t/TQ) = t. Hence, the time t at which the 
behavior switches between approximately adiabatic and 
approximately impulse is ^JtqTq. Once this last "adia- 
batic" instant in the evolution of the system is known, 
interesting features of the post transition state (such as 
the size of the regions ^ in which the order parameter 
is smooth) can be computed. Generalization to other 
universality classes is straightforward and leads to 

r i/{i+uz) ; v/{i+uz) 

where z and v are universal critical exponents (see, e.g.. 

In what follows we adopt this approach to quantum 
systems where dynamics can be approximated as an 
avoided level crossing (see Fig. For simplicity we 
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FIG. 1: Plot (a): structure of energy spectrum (parametrized 
by time) of a generic two- level system under consideration. 
Note the anticrossing at t = and asymptotic form of eigen- 
states: jl), |2). Plot (b): adiabatic-impulse regimes in system 
dynamics. Compare it to a plot of relaxation time scale in 
the Kibble-Zurek theory [fll llH. Note that instants separat- 
ing adiabatic and impulse regimes does not have to be placed 
symmetrically with respect to t = - see Sec. HIT Al and Fig. 

m 

restrict ourselves to two-level systems that possess a sin- 
gle anti-crossing (Fig. ^ in the excitation spectrum and 
a gap preferably going to -t-oo far away from the anti- 
crossing. The latter condition guarantees that asymp- 
totically the system enters adiabatic regime. 

The passage through avoided level crossing can be di- 
vided into an adiabatic and an impulse regime according 
to the size of the gap in comparison with the energy scale 
that characterizes the rate of the imposed changes in sys- 
tem Hamiltonian. If the gap is large enough the system is 
in adiabatic regime, while when the gap is small it under- 
goes impulse time evolution as depicted in Fig. ^b). The 
instants ±£ are supposed to be such that the discrepancy 
between time dependent exact results and those coming 
from splitting of the evolution into only adiabatic and im- 
pulse parts is minimized. Therefore, the AI approxima- 
tion looks like a time dependent variational method where 
i is a variational parameter while adiabatic-impulse as- 
sumptions provide a form of a variational wave function. 

To make sure that assumptions behind the AI ap- 
proach are correctly understood, let's consider time evo- 
lution of the system depicted in Fig. Let the evolution 
start at ti — > —oo from a ground state (GS) and last till 
tf —^ +00. The AI method assumes that the system wave 
function, |^'(i)), satisfies the following three approxima- 
tions coming from passing through first adiabatic, then 
impulse and finally adiabatic regions 

t e [-00, -i] : l^'(i)) w (phase factor)|GS at t) 
t e : l^'(i)) « (phase factor)|GS at - t) 

te[i, +oo] : |(^'(i)|GS at t)p = const. 

Finally, one needs to know how to get the instant t. As 
proposed in [^, the proper generalization of Eq. to 
the quantum case (after rescaling everything to dimen- 
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sionless quantities) reads 



1 



gap(t ) 



at, 



(3) 



where a = 0(1) is a constant. Similarity between Eq. Q 
and Eq. |^ suggests (in accord with physical intuition) 
that the quantum mechanical equivalent of the relaxation 
time scale is an inverse of the gap. 

In this paper we give a simple and systematic way 
for (i) obtaining a analytically; (ii) verification that Eq. 
Q leads to correct results in the lowest nontrivial or- 
der. This method, illustrated on specific examples in Ap- 
pendix^ is based on the observation that a time depen- 
dent Schrodinger equation can be solved exactly in the 
diabatic limit if one looks at the lowest nontrivial terms 
in expressions for excitation amplitudes. This should be 
true even when getting an exact solution turns out to 
be very complicated or even impossible. After obtaining 
a this way, the whole AI approximation is complete in a 
sense that there are no free parameters, so its predictions 
can be rigorously checked by comparison to exact (either 
analytic or numeric) solution. 



III. DYNAMICS IN THE LANDAU-ZENER 
MODEL 

In this section we illustrate the AI approximation by 
considering the Landau-Zener model. In this way we sup- 
plement and extend the results of 0. 

The Landau-Zener model, after rescaling all the quan- 
tities to dimensionless variables, is defined by the Hamil- 
tonian: 




(4) 



where tq is time independent and provides a time scale 
on which the system stays in the neighborhood of an 
anti-crossing. As tq ^ the system undergoes diabatic 
time evolution, while tq 3> 1 means adiabatic evolution. 
There are two eigenstates of this model for any fixed time 
t: the ground state | | (t)) and the excited state | t (0) 
given by: 



T(i)) 
i(i)> 



cos{e{t)/2) sm{e{t)/2) 
-sm{e{t)/2) cos{0{t)/2) 







) 


.12). 



(5) 

where |1) and |2) are time- independent basis states of 
the Hamiltonian Q); cos(6') 



sin(6') 



1/VT 



t 



(6) 



and 6 E [0, tt]. The gap in this model is 
gap 



The instant t is obtained from Equation ijJJ: 
1 

^^^^= = at, 



■i + (i)" 

which leads to the following solution 



t 



1 



(7) 



Dynamics of the system is governed by the Schrodinger 
equation: 

and it will be assumed that evolution happens in the 
interval [ti, tf]. 

As was discussed in fl, by using AI approximation 
one can easily arrive at the following predictions for the 
probability of finding the system in the excited state at 

tf » i 

• time evolution starts at <ti —t from the ground 
state 



(8) 



In this case the system undergoes in the AI scheme 
the adiabatic time evolution from ti to —t, then 
impulse one from — t to t and finally adiabatic from 
£ to > i. 

• time evolution starts at ii = from the ground 
state 



PA7-|(T (<)U(o))P = i(i- 



1 



(9) 



Now evolution is first impulse from t = to t = i 
and then adiabatic from t to tf ^ t. 

In the following we consistently denote excitation prob- 
ability when the system evolves from ti <IC —t, ti = hy 
P, p respectively. Additionally the subscript AI will be 
attached to predictions based on the AI approximation. 
In the first case, the substitution of Q into lO leads 

to 



Pai 



1 



(10) 



Now we have to determine a. It turns out that it can 
be done by looking at the diabatic excitation proba- 
bility. The proper expression and calculation can be 
found in Appendix IXI Substituting rj — 1/2 into (|A6p 
and ljA5|l one gets that to the lowest nontrivial order 
Pai — 1 — t:tq/2, which implies that a — Tr/2 and ver- 
ifies Eq. (O for this case. Note that the calculation 
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FIG. 2; Transition probability when the system starts time 
evolution from a ground state at ti —oo and evolves to 
tf — > +00. Dots: exact expression Ulljl . Solid line: AI predic- 
tion 11211 with Q = 7r/2 determined from diabatic solution in 
Appendix EI Dashed thick line: lowest order diabatic result, 
1 — 'ktq/2, coming from ljA5^ and l|A6|l with r] = 1/2. 



leading to Eqs. (|A5|I and ljA6p is not only much easier 
then determination of exact LZ solution but also 
really elementary. Therefore, we expect that it can be 
done comparably easily for any model of interest. 

Now we are ready to compare our AI approximation 
with a determined as above, to the exact result, i.e., 



P 



exp I 



(11) 



First, the agreement between the exact and AI result is 
up to 0{tq), i.e., one order above the first nontrivial 
term. This is the advantage that the AI approximation 
provides over a simple diabatic approximation performed 
in Appendix ^ Second, we see that the AI expansion 
contains the same powers of tq as the diabatic (small tq) 
expansion of the exact result. Third, Fig. |21 quantifies 
the discrepancies between exact, AI and diabtic results. 
For the AI prediction, we plot in Fig. |21 instead of a 
Taylor series (fTU)! the full expression evaluated in |^ 



Pai 



(arQ)2 + aTQ^ioTQ)^ + 4 + 2 ' 



(12) 



with a = 7r/2. As easily seen the AI approximation 
significantly outperforms a diabatic solution. In other 
words, the combination of AI simplification of dynam- 
ics and diabatic prediction for the purpose of getting the 
constant a leads to fully satisfactory results considering 
simplicity of the whole approach. 

It is instructive to consider now separately three sit- 
uations: (i) dynamics in a nonsymmetric avoided level 
crossing (Sec. IIII All : (ii) dynamics beginning at the anti- 
crossing center (Sec. IIII B|l : and (iii) dynamics starting 
at ti — !■ — oo but ending at the anti-crossing center (Sec. 
IIII C|) . The first case will give us a hint whether 0{t^) 




adiabatic 



FIG. 3: The same as in Fig. 0but for nonsymmetric Landau- 
Zener problem with 5 > 1 - see l|13|l . 



agreement we have seen above is accidental and has some- 
thing to do with the symmetry of the Landau-Zener prob- 
lem. The second problem was preliminarily considered in 
|j9| , but without comparing the AI prediction to exact an- 
alytic one being interesting on its own. Finally, the third 
problem is an example where AI approximation correctly 
suggests at a first sight unexpected symmetry between 
this problem and the one considered in Sec. IIII Bl 



A. Nonsymmetric Landau-Zener problem 

We assume that system Hamiltonian is provided by the 
following expression 



H = 





X = 



1 for t < 
5 for t > 



(13) 



with S > being the asymmetry parameter - see Fig. 
|3fa) for schematic plot of the spectrum. 

Once again, evolution starts at — > — oo from a ground 
state. The exact expression for finding the system in the 
excited eigenstate at the end of time evolution {tf 
+oo) is 



P 



■i7r(l+5)rQ 



■ sinh ( —7^tq5 



T{1/2 + itq6/8) 



Til/2 + ZTQ/8) 



1 r(i - 



,(14) 



and its derivation is presented in Appendix^ Naturally, 
for ^ = 1, i.e., in a symmetric LZ problem, the expression 
(tTIjl reduces to (tTT|l . 

Now we would like to compare (I14|l to predictions com- 
ing from AI approximation. Due to asymmetry of the 
Hamiltonian the systems enters the impulse regime in 
the time interval [—tL,tR] - see Fig. Olb) for illustra- 
tion of these concepts. The instants tj^ and tfi are easily 
found in the same way as in the symmetric case. It is 
a straightforward exercise to verify that according to AI 
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approximation, the probabihty of finding the system in 
the upper state is 



PAi = \{i {-iL)\niR))? 



(15) 



To simphfy comparison between exact (|14|l and ap- 
proximate H15() excitation probabihties we wiU present 
their diabatic Taylor expansions: 



Pai = l-\{l + Vt 



aTQ + 

L(^l + 5){l + V6)\aTQf + 0{T^Q). (16) 

Determination of the constant a is easy and is presented 
in Appendix|X| By putting -q = 1/2 into and HAlOfl . 
one gets that a = n/2. 

The AI prediction can be easily compared to the exact 
result p4|l after series expansion 



64 



-<5)(l + V^)V3 + 0(r|). 



This comparison shows that once again there is perfect 
matching between exact and AI description up to 0{tq), 
i.e., one order better then the simple diabatic approxima- 
tion from Appendix^ despite the fact that we deal now 
with nonsymmetric LZ problem. Moreover, the constant 
a is the same in symmetric and nonsymmetric cases. Fi- 
nally the same powers of tq show up in both exact and 
AI results. 

The nonsymmetric Landau-Zener model is also inter- 
esting in the light of a recent paper [l9| , where it is argued 
that the final state of the system that passed through 
a classical phase transition point, is determined by de- 
tails of dynamics after phase transition. In our system, 
there are two different quench rates: tq before the tran- 
sition and Tq = tq6 after the transition. Substitution 
oi 5 = Tq/tq into (|14|) shows that the final state of non- 
symmetric Landau-Zener model expressed in terms of the 
excitation probability depends on both tq and Tq , so that 
it behaves differently than the system described in . 



B. Landau-Zener problem when time evolution 
starts at the anti-crossing center 

Now we consider the LZ problem characterized by 
Hamiltonian in the case when time evolution starts 
from the ground state at the anticrossing center, i.e., 
t, = 0. It means that |5'(0)) cx (|1) - |2)) /V2. Exci- 
tation probability at — > -l-oo equals exactly (see Ap- 
pendix 



sinh ( — —\ e 

tttq V 4 / 



r 1 



8 / 




(17) 



From the point of view of AI approximation the evo- 
lution is now simplified by assuming that it is impulse 
from ti = Q to t and then adiabatic from t to tf ^ +oo. 
This case was discussed in 'o^, and it was shown that 
substitution of Q into © leads to a simple expression 



_ 1 1 / 2 

" 2 " 2 Y ^ " (arQ)2 + aTg^iaTg)^ + 4 +~2' 



expanding it into diabatic series one gets 



(18) 



PAi = \- ^V^+lic^TQ^/' + 0{t'J'). (19) 



Determination of a constant a is straightforward (see Ap- 
pendix 0. Putting 77 = 1/2 into (|A8IA9|I one easily gets 
a = 7r/4 and confirms that the first nontrivial term in 
(|19|l is indeed oc ^/tq- This value of a (« 0.785) is in 
agreement with the numerical fit done in |^ where opti- 
mal a was found to be equal to 0.77. Small disagreement 
comes from the fact that now we use just the ^ — ^^arg 
part for getting a, which is equivalent to making the fit 
to exact results in the limit of tq ^ 1. In 9] the fit of 
the whole expression (|18|) in the range of < rg < 6 was 
done. 

Having the exact solution Ijl?!) at hand, we can verify 
rigorously accuracy of H18I19II with the choice of a = 7r/4. 
Expanding (|17|) into diabatic series one gets 



1 



^3/2 / 4ln2^ 
64 



o{4 



(20) 

As expected the first nontrivial term is in perfect agree- 
ment with H19|l once a = n/A. On the other hand, the 

term proportional to Tq^ is about 10% off. Indeed, af- 
ter substitution a = 7r/4 into ifTOI) one gets the third 
term equal to tt^^'^Tq^ /64:, which has to be compared to 

w IAtt^/'^Tq^'^ /64: from It shows that the AI ap- 

proximation in this case does not predict exactly higher 
nontrivial terms in the diabatic expansion. 

Nonetheless, expression (|18|l works much better then 
the above comparison suggests. It not only signifi- 
cantly outperforms the lowest order diabatic expansion, 
i — J -y/TTTg , but also beautifully fits to the exact result 
pretty far from the rg ^ 1 regime. Fig. 0] leaves no 
doubt about these statements. It suggests that AI pre- 
dictions should be compared to exact results not only 
term by term as in the diabatic expansion, but also in a 
full form. 



C. One half of Landau-Zener evolution 

In this section we consider exactly one half of the 
LZ problem. Namely, we take the Hamiltonian Q and 
evolve the system to the anti-crossing center, tf = 0, 
while starting from the ground state at ti —f —00. 
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FIG. 4: Transition probability when the system starts time 
evolution at the anti-crossing center from a ground state and 
evolves to tf ^ +oo. Dots: exact expression H17|l . Solid 
line: AI prediction 1181 with a = tt/4 determined from the 
diabatic solution in Appendix^ Dashed thick line: lowest 
order diabatic prediction 1/2 — ^-ktq /4 coming from IIA8IA9II 
with 77 = 1/2. 



Let's see what the AI approximation predicts for exci- 
tation probability of the system at = 0. According to 
AI, evolution is adiabatic in the interval [—00,—t] and 
then impulse in time range [— 0]. It implies that 

|(t (o)|*(o))p « |(T m^i-i))\' « KT (0)1 i 



A simple calculation then shows that 



(21) 



A quick look at (O shows that the excitation probability 
for one-half of the Landau-Zener problem is supposed to 
be, according to the AI scheme, equal to the excitation 
probability when the system starts from the anti-crossing 
and evolves toward tf +00. 

To check that prediction we have solved the one-half 
LZ model exactly, see Appendix ^ and found that in- 
deed both probabilities are exactly the same, so the AI 
approximation provides us here with a prediction that is 
correct not only qualitatively but also quantitatively. 



IV. DYNAMICS IN A GENERAL CLASS OF 
QUANTUM TWO LEVEL SYSTEMS 

To test predictions coming from the AI approximation 
on two level systems different than the classic Landau- 
Zener model, we consider in this section dynamics in- 
duced by the Hamiltonian 



H 



sgn(i) 



1 



-sgn(i) 



(22) 



where ij G (0, 1) is a constant parameter (when 77 = 1/2 
the system reduces to the Landau-Zener model). The 
eigenstates of Hamiltonian l|22|l are expressed by formula 
© with 



cos(6') = 



sgn(t) 



and 



sin(6») 



2t? 
1-17 



Extending the notation of Sec. Illll to the system (|22|l one 
has to replace definition © by 



e = sgn(i) 



t 



and then expressions (|8I9|I providing excitation probabil- 
ities are valid also for the model (|22|l . 

The instant t for the Hamiltonian l|22|l is determined 
by the following version of Eq. ^ 



1 



2ri 



at. 



(23) 



It seems to be impossible to obtain solution of (|23|l ex- 
actly. Fortunately, a diabatic perturbative expansion for 
t can be found - see Appendix lO Once the diabatic ex- 
pansion of t is known, we can present different predictions 
about dynamics of the model H22|) . 

First, we consider transitions starting from a ground 
state at ti —00 and lasting till tf +00. Using the 
results of Appendix [CI one gets the following expression 
for the excitation probability: 



Pai 



1 + 



l-(aTQ)2" + (l-,y)(«rQ) 



4r; 



i(l - 77)(2 - 3r;)(aTQ)6" + 1(1 - ^)(1 - 2r;) 



{3-Arj){arQr + 0{r'J") 



(24) 



This result is equivalent to (|10|) when one puts rj = 1/2. 
The constant a was found in Appendix ^ from the dia- 
batic solution and equals IIA6|l 



a = (1 — ?7)r(l — ri)i . 

The exponent 2ri in the lowest order term in H24|l was pos- 
itively verified ( Appendix 0. Additionally, we have per- 
formed numerical simulations for 77 — 1/3,2/3,3/4,5/6. 
A fit to numerics in the range of tq <C 1 has confirmed 
both the exponent 2ri and the value of a. The comparison 
between the AI prediction for rj = 1/3, 2/3 and numerics 
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is presented in Fig. |Sfa) . The lack of exact solution for 
this problem does not allow us for a more systematic an- 
alytic investigation of the AI approximation in this case. 
Nonetheless, a good agreement between AI and numerics 
is easily noticed. Once again, the AI prediction outper- 
forms the lowest order diabatic expansion. 

As before, we would like to consider transitions starting 
from the anticrossing center, i.e., ti = 0. In this case the 
excitation probability equals 



PAI 



VT 



^-^(«TQr + i(i-7?) 



O 



(25) 



As above, the diabatic calculation of Appendix ^ veri- 
fies the exponent of the first nontrivial term in 1251) and 
provides us with the following prediction for a (|A9|I 



{1 - ri)T {1 -r])'" cos (^^riY . 



After setting = 1/2 formula (|25|l becomes the same as 
(|19|l . Due to the lack of exact results we have carried out 
numerical simulations for rj — 1/3,2/3,3/4,5/6. As ex- 
pected, this calculation has confirmed that the exponent 
of tq is indeed ij for small tq and that a is provided by 
()A9|I . The AI approximation vs. numerics and diabatic 
expansion for 77 ~ 1/3,2/3 is presented in Fig. EJb). 
Once again an overall agreement is strikingly good. 



V. QUANTUM PHASE TRANSITION IN ISING 
MODEL VS. ADIABATIC-IMPULSE APPROACH 

In this section we will illustrate how ideas coming from 
the adiabatic-impulse approach can be used in studies of 
dynamics of quantum phase transition. As an example 
we choose the quantum Ising model recently considered 
in this context H, Q. The quantum Ising model is 
defined by the following Hamiltonian: 



N 



'n'-'n+lj I 



(26) 



where N is the number of spins. The quantum phase 
transition in this model is driven by the change of a di- 
mensionless coupling g. The transition point between 
paramagnetic phase {g > 1) and ferromagnetic one (0 < 
5 < 1) is at (7 = 1. To study dynamics of (|26|l one as- 
sumes that the system evolves from time i' = — 00 to 
time t' = 0, and takes 
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where Tq provides a quench rate. The quantity of interest 
is density of topological defects (kinks) after completion 




FIG. 5: Numerics: bold dashed line, AI approximation: solid 
line, lowest order exact diabatic expansion f Appendix IXl: 
dotted line, (a): evolution starts from ti — > —00. (b): evolu- 
tion starts from ti = 0. In both (a) and (b) plots evolution 
ends in tf ^ +00. 



of the transition, i.e., at 5 = 0. It equals 



N-l 



= lim / ^— (1 — cr^cr^ I 1 ) 



(27) 



The quantum Ising model obviously possesses 2^ dif- 
ferent energy eigenstates so it seems to be hopeless to 
expect that the two level approximation would be suffi- 
cient. Therefore, it is a remarkable result of Dziarmaga 
j3| , that dynamics in this system can be exactly described 
by a series of uncoupled Landau-Zener systems. Due to 
lack of space, we refer the reader to 0, and present just 
the main results and their AI equivalents. 

The density of defects (|27|> in Dziarmaga's notation 
reads as 

27r . _^ 
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where pk 's are defined as 

Pk = I cosik/2)ukit' = 0) - sin(fc/2)wfe(i' = 0)p, (28) 
with Uki Vk satisfying the following Landau-Zener system 

(29) 




t - 4 (t' + 



> fc) sin k, 



that \uf{tf)\ « |u/(+oo)|. Initial conditions (|30|l mean 
that the system starts time evolution at t = — c» from ex- 
cited state and we are interested in probability of finding 
it in the ground state at i/. Elementary algebra based 
on AI approximation shows that this probability equals 

K(i/)P = |(T {-i)\i{i))? = \{i M)IT(<))P. 

Therefore, it is provided by (|HJl and (|12|l . For any fixed e 
consistent with lowest order approximation of sin k and 
cos k one gets 



'^Q = 4tq sin^ fc, 

where changes from —oo to 0. The initial conditions 
for that LZ system are the following 



0.172- 



|ufc(<' = -oo)| = 1 , Vk{t' = -oo) = Q. 



(30) 



Due to the symmetry of the whole problem, it is possible 
to show that pk = p-k, and therefore it is sufficient to 
restrict to k > from now on. 

Evolution in H29|l lasts from ti ~ — oo to tf — 
2TgSin(2/c). Since obviously tf does not go to +oo, it is 
important to determine where i/ is placed with respect 
to ±i. 

We solve the equation t = \tf\ (where we substitute 
k = kc). After using Q with a = 7r/2, and simple 
algebra one arrives at 



1 



- 1 = 



cos kr 



sin kr 



which leads to the result 



sin kr — 1 — 



1 



Defining kc to be between and 7r/2 and doing some 
additional easy calculations one gets that when t'q > ^ 
we have: (i) tf > i for k £ [0,kc]; (ii) — < < i/ < t for 
k G {kc,TT — kc); (iii) tf < — < for A; e [tt — fcc,7r]. Moreover, 
when Tq < ^ one can show that \tf\ < i for any fc's of 
interest. It means that the evolution ends in the impulse 
regime for small enough Tq. 

Now, as in 0, IE 01 J we consider adiabatic time evolu- 
tions of the Ising model. In our scheme it clearly corre- 
sponds to Tq 3> ^ . As the system undergoes slow evolu- 
tion it is safe to assume that only long wavelength modes 
are excited, which means that fc <C 7r/4 are of interest. 
This allows to approximate sin fc ~ fc and cos fc ~ 1 , 
which implies that H28|l turns into pk ~ l^^fc(^/)P and 
that 

1 f 

n« - / dk\ukitf)\^, 

Jo 

with e <C 7r/4. Since fc g [0,e] corresponds to the above 
mentioned (i) case, tf > t , the AI approximation says 



(31) 



where the prefactor, was found from expansion of the 
integral into 1/wr^ adiabatic series with +oo at 



adiabatic series with Tq 

fixed e. In the derivation of this result sin^ fc in the ex- 
pression for Tq was approximated by fc^. 

First of all, the prediction H31() provides correct scal- 
ing of defect density with Tq. The prefactor, 0.172, has 

to be compared to ^ « 0.16 (exact result from Q and 
numerical estimation from f^l) or 0.18 (approximate re- 
sult from Q). Our prefactor matches these results re- 
markably closely concerning simplicity of the whole AI 
approximation. Slight overestimation of the prefactor in 
comparison to exact result comes from the fact that AI 
transition probability, Eq. H12I) . overestimates the exact 
result, Eq. l(TT|) . for large tq's (Fig. It is interesting, 
to note that in the Kibble-Zurek scheme used for de- 
scription of classical phase transitions an overestimation 
of defect density is usually of the order of 0(1) 2. 12 , 13j |, 
while here discrepancy is smaller than 10%. 

Therefore, the AI approximation provides us with two 
predictions concerning dynamics of the quantum Ising 
model: (i) it estimates when evolution ends in the asymp- 
totic limit; (ii) it correctly predicts scaling exponent and 
number of defects produced during adiabatic transitions. 
The part (ii) can be calculated exactly in the quantum 
Ising model H26|l . Nonetheless, in other systems undergo- 
ing quantum phase transition the two level simplification 
might be too difficult for exact analytic treatment, e.g., 
as in the system governed by (|22|l . Then the AI analy- 
sis might be the only analytic approach working beyond 
the lowest order diabatic expansion. Besides that the AI 
approach provides us with a quite intuitive description 
of system dynamics, which is of interest in its own. Es- 
pecially, when one looks at connections between classical 
and quantum phase transitions. 



VI. SUMMARY 

We have shown that the adiabatic-impulse approxima- 
tion, based on the ideas underlying Kibble-Zurek mecha- 
nism |l(]llll| . provides good quantitative predictions con- 
cerning diabatic dynamics of two level Landau-Zener like 
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systems. After supplementing the splitting of the evo- 
lution into adiabatic and impulse regimes by the exact 
lowest order diabatic calculation, the whole approach is 
complete and can be in principle applied to different sys- 
tems possessing anti-crossings, e.g., those lacking exact 
analytic solution as the model H22|l . We expect that the 
AI approach will provide the link between dynamics of 
classical and quantum phase transition thanks to usage 
of the same terminology and similar assumptions. 

We also expect that different variations of the classic 
Landau-Zener system discussed in this paper can be ex- 
perimentally realized in the setting where the sweep rate 
can be manipulated. Indeed, the model H22() arises once 
a proper nonlinear change of external system parame- 
ter is performed. The experimental access to studies of 
nonsymmetric Landau-Zener model can be obtained by 
change of the sweep rate after passing an anti-crossing 
center. The case of evolution starting from or ending at 
the anticrossing center can also be subjected to experi- 
mental investigations, e.g., in a beautiful system consist- 
ing of the smallest available quantum magnets (cold Fcg 
clusters) 0- 

We are grateful to Uwe Dorner for collaboration 
on nonsymmetric Landau-Zener problem, and Jacek 
Dziarmaga for comments and useful suggestions on the 
manuscript. Work supported by the U.S. Department 
of Energy, National Security Agency, and the ESF 
COSLAB program. 



to +00. The initial condition is such that Ci(— 00) = 1 
and C2{—oo) = 0. The simplification comes when one 
assumes a very fast transition, i.e., tq — > 0. Then it is 
clear that Ci{t) = 1 + 0{t^) with some /3 > 0. Putting 
such Ci(i) into (|A3|) one gets 



C2(+oo) = — J da;exp(-i(l-77)|a;|i-'-) -h 

(higher order terms in tq), (A4) 

which after some algebra results in 

P=|Ci(+oo)|2 ^ 1 - \C2{+oo)\^ ^ 1 - {arg f^^ + 

(higher order terms in tq), (A5) 



where 



a= {1 - r])T {1 - r])~ 



(A6) 



Time evolution starting from the ground state 
at ti = (anti-crossing center): since initial wave 
function is (|1) - |2))/\/2 we have Ci(0) = 1/V2 and 
^2(0) — —I/V2 and we evolve the system till tj ^ +00 
with the Hamiltonian H22|) . For fast transitions one has 
that Ci{t) = 1/^/2 + O(t^), where /3 > is some con- 
stant. Integrating (jA3|) from to -foo one gets: 



APPENDIX A: EXACT DIABATIC 
EXPRESSIONS FOR TRANSITION 
PROBABILITIES 



C2(+00) + i= 



'Q 



+00 



dxexp I — i(l — ri)x'^ 



2V2i JO 

-|- (higher order terms in tq), (A7) 



We would like to provide lowest order exact expressions 
for transition probabilities in a class of two level systems 
described by the Hamiltonian H22() . 

First, we express the wave function as 



\^{t)) = Ci(i)exp 



-(^-j^l^ ID 



C2(i)exp|^^^^fc^||2), (Al) 

were exponentials are =Fz f dtsgn(i)|t/rQ|T~^/2. Within 
this representation dynamics governed by the Hamilto- 
nian (|^ reduces to 

.6,,), ftl2exp( -"-i'l'^ | (A2) 



which can be easily shown to lead to 



1 1 



|Ci(-foo)|^ = 1 - |C2(+(^)|^ = - - -(arg)" -I- 



(higher order terms in tq), 



where 



(1 - ri)T {I -Tj)^ cos (^'^rij" . 



(A8) 



(A9) 



Time evolution from 00 to tf +00 in 

the nonsymmetric LZ model: we assume that the 
Hamiltonian is given by (|22l for i < 0, while for t > it 
is given by 1)22(1 with tq exchanged by tqS where 6 > 
is the asymmetry constant - see ((13|l for 77 = 1/2 case. 
Integration of l|A3p separately in the intervals [— oo,0] 
and [0, -t-00] leads to the following prediction 



■r u\ / -»(l-77)|t|— , 
«<^2(i) = exp ^ 3- I (A3) 

Time evolution starting from the ground state 

at — > — 00: we would like to integrate l|A3|l from —00 



P = |Ci(+oo)|2 ^ 1- |C2(-f^)|' 

(1 + 5")' (arQ)2" -f 
(higher order terms in tq). 



= 1- - 
4 



with a given by (|A6I 



(AlO) 
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APPENDIX B: EXACT SOLUTIONS OF THE 
LANDAU-ZENER SYSTEM 

In this Appendix we present derivation of exact solu- 
tions of the Landau-Zener problem that are used in the 
main part of the paper. The general exact solution for 
Landau-Zener model evol ving according to Hamiltonian 
Q was first discussed in uM and then elaborated in a 
number of papers, e.g., |20|. Here we will use that gen- 
eral solution for getting predictions about different time 
evolutions. The problem is simplified by writing the wave 
function as (|A1(I with ?/ = 1/2. Then one arrives at the 
equations l|A2|l and ljA3|l once again with rj = 1/2. Com- 
bining the latter ones with the substitution 



one gets 



where 



U2{Z) 



U2{t) = C2(t)e^-« 



fc - ^ + ^ ) U2{Z) - 0, 



4 



(Bl) 



The general solution of this equation is expressed in terms 
of linearly independent Weber functions D^k-iiHz) 0, 
l2l] |. By combining this observation with rj = 1/2 version 
of HA3|I one gets 



[aD^k~i{iz) + bD. 



k-iy 



+ [aD-k-i{iz) + hD-k-i{~iz)]\2), 



-^z)]|l) 
(B2) 



where z is defined in I|B1(I . Though there is a simple one- 
to-one correspondence between z and we will use both 
z and t in different expressions to shorten notation. 

To determine constants a and b in different cases one 
has to know the following properties of Weber functions 
[2l| . First, V \arg{s)\ < 37r/4 one has 

Dr^is) - e-'^'/V"[l + 0(s-2)]. (B3) 

Second, V— 57r/4 < arg{s) < — 7r/4 

2^ 



An(s) = e' 
Third, as s 



one has 



r(— to) 



r(l/2- to/2) 



g-i(m+l)7r/2^_ 



Ois) 



Finally, 



^Drn{s) = mDm-l{s) - \sDm{s). 

as 2 



i{is). 
(B4) 



(B5) 



(B6) 



Exact solution of nonsymmetric LZ problem: 

the evolution starts at ti — s- —oo, i.e., up to a phase 
factor |\['(— oo)) ~ |1). The Hamiltonian is given by 
(P)l . Using P3IB4IB6P one finds that a = and 
b — y^exp(— 7rTQ/16)/2. Substituting them into HB2|) 
results in 



g-'""Q/16gi37r/4 

[(fc + l)D-k-2Hz) - tzD-k-i{-tz)]\l) 



\^it<0)) 



+^e--«/i6i?_,_i(-*z)]|2). (B7) 
Using IjBSp one finds that this solution at t = becomes 



l*(0)) 



— p-TTTQ/ie i37r/4_ 



V^2-'=/2 



r(i/2 + fc/2)' ' 

2 ^ y2T{l + k/2y'- 



For t > the Hamiltonian changes its form, see (|13|) , and 
one has to match (jB8|l with (|B2|I having tq replaced by 
tqS. Substitution of 



-tttq/IG 



T{l + kS/2) jl_ 
T{l + k/2) \d 



T{l/2 + kd/2) 

' r(i/2 + fc/2) 



4 ' 
r(l/2-| 



-fc(l-<5)/2 



kS/2) 



T{l + k6/2) 
r(l + k/2) 



r(l/2 + fc/2) 



into l|B2p gives the wave function > 0)). Taking 
one minus squared modulus of the amplitude of finding 
the system in the state |2) one obtains the excitation 
probability of the system at tf +oo in the form (|14|) . 

Exact solution of LZ problem when evolution 
starts from a ground state at anticrossing center: 
the Hamiltonian of the system is given by Q and we 
look for a solution that starts at t = from the (ground) 
state: |^'(0)) = (|1> - |2))/V2. The constants a and b 
from HB2II turn out to be equal to 



^fc/2 r(i + fc/2) 

20F ' 



2fe/2r(l/2 + A:/2)e"/4 
4^2^ 



2'-/^r(l/2 + fc/2)eW4^ ^ 2,/2r(l+_fc/2) 



4V27r 



2V^ 



Having this at hand, it is straightforward to show that 
one minus squared modulus of the amplitude of finding 
the system in the state |2) at tf — s- -|-oo, i.e., excitation 
probability, equals ((T7|l . 
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■^Ul) ^L(2r'^U(2)^U(l) 

FIG. 6; Schematic plot of the recurrence method for getting 
solution of IIC2II . Solid lines x ^—n and (i — x, dashed line is 
a construction of a recurrence solution (O and (EH. The 
plots are in a /3 x /3 box. 



One-half of the Landau-Zener problem: the ex- 
citation probability of the system that started time evo- 
lution at ti = —oo from ground state, and stopped evolu- 
tion at t/ = 0, turns out to be equal to l|17|) . The simplest 
way to prove it comes from an observation that this prob- 
ability equals one minus the squared overlap between the 
state (|1) - |2))/\/2 and ||B8)) . Then, straightforward cal- 
culation leads immediately to expression H17|) . 



APPENDIX C: SOLUTION OF EQ. 

Equation (|23|l can be solved in the following recursive 
way. First, one introduces new variables: 



1 



(CI) 



In these variables Equation H23() becomes: 

a;T^=/3-a;, (C2) 

and we assume that /3 > 1 (diabatic evolutions). 

A quick look at Fig. El makes clear that the solution 
can be obtained by considering a series of inequalities, 
numbered by index i, in the form 

XL(i) <X< Xu(i), (C3) 

where both lower {x^^i)) and upper (x^(i)) bounds satisfy 
the the same recurrence equation 



f3-[l3- XL^ua^i) 



1-') 



(C4) 



with the initial conditions 0:^(0) = Oi 2:^;7(o) = P- Consid- 
ering a few iterations one can show that the solution of 
T2|) can be conveniently written as 



1 



X 1-1 



[l-{arf^' + {l-7j){aTQ)^^+ 



i(3ry - 2)(1 - vK^tq)^'' + ^(1 - r;)(l - 2r,) 
(3-4,?)(arQ)«'' + 0(r^"")]. (C5) 



Using (|C5I) . different quantities of interest can be deter- 
mined, e.g.. 



1 



2r, 



l-,7(aTQ)2'' + i,;(l-^)(«rQ)4"- 



t = 



(l-377)(2-3r;)(arQ)«''-|-0(r^"'') 



(C6) 
(C7) 



Substitution of r/ = 1/2 into (|C6p and (|C7|) reproduces di- 
abatic series expansions of i, t from the standard Landau- 
Zener theory iQ. 
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